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The storage properties of methane gas in Vycor porous glass (5.7 nm) are characterized in a wide 
pressure range from 0.7 MPa-89.7 MPa using Nuclear Magnetic Resonance (NMR). We demonstrate 
the capability of high held NMR relaxometry for the determination of the methane gas storage 
capacity and the measurement of the Hydrogen Index, to a high degree of accuracy. This helps 
determine the excess gas in the pore space which can be identified to exhibit Langmuir properties in 
the low pressure regime of 0.7 MPa to 39.6 Mpa. The Langmuir model enables us to determine the 
equilibrium density of the monolayer of adsorbed gas to be 8.5% lower than that of liquid methane. 

We also identify the signatures of multilayer adsorption at the high pressure regime from 39.6 Mpa to 
89.7 Mpa and use the Brunauer-Emmet-Tellcr (BET) theory to determine the number of adsorbed 
layers of methane gas. We show how these measurements help us differentiate the gas stored in 
the Vycor pore space into free and adsorbed fractions for the entire pressure range paving way for 
similar applications such as studying natural gas storage in gas shale rock or hydrogen storage in 
carbon nanotubes. 


I. INTRODUCTION 

Nuclear magnetic resonance (NMR) is a well- 
established technique for probing the dynamics of flu¬ 
ids in porous media with disordered geometries, such 
as rocks and biological systems [1, 2]. NMR investiga¬ 
tions of fluids in porous media are generally performed 
through measurements of relaxation times and diffusion 
coefficients. This enables the determination of the fluid 
properties such as viscosity, chemical composition and 
fluid-pore surface interactions. In addition, the geomet¬ 
rical properties of the porous matrix such as tortuosity 
(a), surface to volume ratio (S/V) and diffusive perme¬ 
ability may also be determined by studying the relaxation 
and time-dependent diffusion coefficient of the confined 
fluid [3-5]. For these reasons, NMR has become an in¬ 
valuable tool, especially in oil and gas exploration appli¬ 
cations where it is routinely applied for downhole logging 
and laboratory rock core analysis. NMR has also found 
broader applications in energy research including investi¬ 
gations of methane gas in shale rocks [6-8] and hydrogen 
storage in carbon nanotubes and other nanoporous ma¬ 
terials [9, 10]. 

The storage of gas molecules in porous media has seen 
renewed interest especially due to recent advances in nat¬ 
ural gas production from gas shale reservoirs. Gas shales 
are organic-rich mudstones characterized by the presence 
of nanometer size organic kerogen pores, hosting natural 
gas. The gas molecules in the nanopores coexist in two 
phases: free gas that has bulk-like properties and ad¬ 
sorbed gas that undergoes strong interactions with the 
pore walls [8] . These two phases are in fast exchange dur¬ 
ing the time scales of the NMR experiment [11]. In addi¬ 
tion, the nature and dynamics of these two phases depend 


on the temperature and pressure of the system, making 
adsorption studies of great importance for understand¬ 
ing natural gas storage and extraction. Though NMR of 
fluids in porous media is a well-established field covering 
studies such as those of diffusion of supercritical fluids in 
nanoporous media using pulsed field gradient techniques 
[12-14] and methane gas adsorption on graphite [15], the 
storage properties of natural gas in nano pores at pres¬ 
sures close to the highest pressures encountered in gas 
shale reservoirs (~ 90 MPa) have not been studied to 
our knowledge. 

In this work we describe a detailed study of methane 
gas storage in nanoporous Vycor glass. We confirm spin- 
rotation to be the dominant contributor to the relax¬ 
ation Hamiltonian of bulk methane gas, up to ultra-high 
pressures of 70 MPa. We demonstrate the capability of 
high field NMR to determine the gas capacity, Hydrogen- 
Index (HI) (sec. Ill C) , adsorbed gas density, and enable 
the separation of free and adsorbed gas quantities in Vy¬ 
cor nanoporous glass at a wide range of pressures (0.7 
MPa-89.6 MPa). We show how this analysis is enabled 
by the unique spin-spin relaxation times of the methane 
gas in the pore space. The quantification of methane gas 
stored in the nanopores together with the pore volume 
allow for a direct measurement of the effective Hydrogen- 
Index. We determine the excess gas in the pore space due 
to adsorption by using either monolayer or multilayer ad¬ 
sorption models for the appropriate low pressure and high 
pressure regimes. This enables the determination of the 
equilibrium density of the adsorbed layer and the number 
of adsorbed layers permitting the separation of the total 
gas into free and adsorbed fractions. 
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FIG. 2. Two dimensional T\ — T 2 correlation pulse sequence 
used in this work. As discussed in the text, the time interval 
r in the indirect dimension was varied from 10ms to 60s. The 
number of n pulses was 32,000 and the time echo time T was 
200/rs. 


FIG. 1. Schematic representation of the experimental setup 
used in this work. The apparatus is composed of two pressure 
generators, the hand crank pump and the syringe pump. The 
upper pressure limit for this setup is 206.8 MPa. The pressure 
is regulated digitally to within 0.35 MPa. 

II. EXPERIMENTAL METHODS 

All experiments were performed on a Bruker 400MHz 
system equipped with high-pressure apparatus as shown 
in the schematic in Fig. 1. The main pressure genera¬ 
tion was carried out through the syringe pump and the 
hand crank pressure generator, which were connected to 
the high pressure NMR sample tube. The hand crank 
pressure generator had a volume of 18 ml and was pur¬ 
chased from High Pressure Equipment Company (HIP). 
The syringe pump (Daedalus Innovations LLC) with a 
capacity of 6.7 ml and pressure rating of 413.7 MPa was 
used for secondary compression. The NMR sample tube 
(Daedalus Innovations) made of zirconium, had outer and 
inner diameters of 5.0 mm and 3.0 mm respectively, and 
was 92 mm long. Methane gas of 99.9% purity was pur¬ 
chased from Air-Gas (Radnor Township-PA) in a cylinder 
pressurized to 13.8 MPa. A three way valve was used to 
direct the methane gas to either the sample tube or the 
vacuum pump. All the connections were well evacuated 
before the experiments to purge out all oxygen in the 
system and avoid any additional relaxation effects. The 
high pressure components were connected with pressure 
tubing and valves rated to 413.7 MPa. 

Vycor nanoporous glass manufactured by Corning Inc 
(New York, USA), was purchased in the form of cylin¬ 
drical rods of diameter 2.96 mm. They were cut into 
samples of 11.50 mm to fit in the center of the NMR 
coil (16.68 mm). The pore size of the Vycor rods was 
5.7 nm and the porosity 28% of the sample volume, as 
determined by BJH adsorption studies. 

The transverse relaxation times (T 2 ) were measured 
using the Carr-Purcell-Meiboom-Gill (CPMG) pulse se¬ 
quence [16, 17]. The two dimensional T 1 -T 2 correlation 


experiments were performed by using the inversion re¬ 
covery pulse sequence in the indirect dimension and the 
CPMG pulse sequence in the direct dimension, as shown 
in Fig. 2. The time interval t for encoding in the in¬ 
direct dimension was varied from 10ms to 60s (the time 
at which the magnetization is fully relaxed). The num¬ 
ber of 7T pulses was chosen to be 32,000 to cover the en¬ 
tire T 2 spectrum and the echo time T was 200/rs. Diffu¬ 
sion effects were confirmed to be negligible at such short 
echo times by observing the attenuation of the first five 
echoes at various echo times [18]. The last delay (repe¬ 
tition time) between scans was varied for each pressure 
from 20s to 60s. The temperature was kept constant at 
22° C throughout all the experiments with an accuracy of 
±1°C. 

The RF probe and the receiver were calibrated to mea¬ 
sure the number of l H nuclei in the sample by carrying 
out T 2 measurements of 5/iL to 40/ii of H 2 0 samples at 
22°C and atmospheric pressure. We determined the ac¬ 
curacy of the calibration to be within 4% by comparing 
the NMR signals of bulk methane gas at various pres¬ 
sures with its density values from National Institute of 
Standards and Technology (NIST) [19]. 

III. RESULTS AND DISCUSSION 

A. Relaxation of free methane gas 

The transverse ( T 2 ) and longitudinal (T)) relaxation 
times of bulk methane gas were measured for a wide 
range of pressures from 3.5 MPa to 68.9 MPa at 22°C. 
The critical pressure and temperature of methane are 4.6 
MPa and -82.75°C respectively, and therefore most ex¬ 
periments were conducted in the super critical state. The 
measured T\ relaxation times are consistent with the lit¬ 
erature data at the low pressure range measured by Wat¬ 
son et al. [20] as shown in Fig. 3. The relaxation times 
of bulk methane gas are dominated by the spin rotation 
mechanism at low to moderate pressures while at high 
pressures, due to the increased density, the mean free 
path decreases and inter molecular dipolar interactions 
could start playing a role [21]. 
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FIG. 3. Longitudinal (Ti) and transverse (T 2 ) relaxation 
times of bulk methane gas at various pressures at 22° C. To 
convert the magnetization decays into relaxation distribu¬ 
tions, a fast inverse Laplace transform algorithm was used 
[22] and each experimental point corresponds to the log mean 
of the entire ID distribution. The reduced value of T 2 in com¬ 
parison to Ti is because of the interactions with the sample 
tube walls. The literature values of Tuiterature (magenta) are 
taken from Ref. [20]. 


In the motional narrowing limit the relationship be¬ 
tween the relaxation rates is given by, 


! , 2 (Cm + 2C\)I\Tf 

Tr 1 = t. r 1 = —— kT . (i) 


In the above expression Tp is the rotational correlation 
time, k is the Boltzmann’s constant, I\ is the moment of 
inertia of the spherical molecule, T is the temperature 
and C ||, C± are the principal components of the spin 
rotation tensor [23]. The correlation time, Tp is given by, 


3/i D 

TF ~ Aa 2 kT' 


( 2 ) 


where D denotes the diffusion coefficient of the molecule 
which depends on the viscosity as D = kT/Omar], with a 
being the radius of the molecule and 7 the viscosity. 

Bulk methane gas is in the motional narrowing limit at 
the Larmor frequency of the experiments (/ = 400 MHz) 
and therefore T\ should be equal to T 2 . However, as 
shown in Fig. 3 the T 2 values are consistently shorter 
than the T\ values at all measured pressures. This is 
because the spin-spin relaxation times (T 2 ) are sensitive 
to the slow motions of methane molecules caused by their 
interaction with the NMR sample tube walls. This may 
be understood by inspecting the frequency dependence 
of the relaxation rates [24], 



37WW±1) |j(1)h 


+J {2) ( 2w)], 


( 3 ) 


and, 

(4) 

where p 0 is the vacuum permeability, I is the spin 
number, 7 is the gyromagnetic ratio and r is the inter- 
nuclear vector between the 1 H nuclei. The spectral den¬ 
sities J (w) can be obtained by the Fourier transform 
of the autocorrelation function Q{t) = (. B(t)B(t + r)) 
which describes the time dependent fluctuations of the 
local magnetic field B(t). The T 2 relaxation times are 
dominated by the = 0 ) term and are therefore 

very sensitive to the low frequency or slow motions while 
the Ti relaxation times are sensitive to the much higher 
Larmor frequency (ujq & 2w 0 ). 

Assuming Gaussian diffusion in ID, the mean displace¬ 
ment of the methane gas molecules may be calculated as, 
(A x(t)) = (2DI) 1 / 2 , where D is the bulk diffusion coeffi¬ 
cient of methane and t is the diffusion time period. At a 
pressure of 20.7 MPa the spins have a self-diffusion coef¬ 
ficient of 7.7 • 10~ 8 m 2 /s [25] and a spin lattice relaxation 
time of 5.2s [26], resulting in displacements of 0.9 mm 
during their T) lifetimes. As the high pressure tube has 
an inner diameter of 3.0 mm, a fraction of the spins inter¬ 
act with the tube walls. This makes the magnetization 
decay of the T 2 measurement (CPMG sequence) multi¬ 
exponential. The T 2 relaxation times shown in Fig. 3 are 
the log mean values of the T 2 distributions. The spin lat¬ 
tice relaxation times (Tf), being only sensitive to motions 
that occur around the Larmor frequency (/ = 400 MHz), 
are not affected and the corresponding magnetization de¬ 
cays are mono exponential. 

The total Hamiltonian for the spin ensemble in the ro¬ 
tating frame may be written as, TLtotai = Hfj 0 + 71 sr 
where T~Lsr describes the spin rotation interaction and 
T~Lfj D is the dipolar Hamiltonian respectively. At the 
pressure of 68.9 MPa the density of bulk methane is still 
26% lower than that of liquid methane. At low pressures 
where the total Hamiltonian is completely dominated by 
spin rotation interaction an empirical relation between 
the spin lattice relaxation times and the methane gas 
density was obtained [28], 


Tisr — c ■ p ■ T -3 / 2 , (5) 

where p is the density in g/cm 3 , T is the temperature 
in K and the constant c = 1.57• 10 5 [cm 3 /\ 3 / 2 sg -1 ]. This 
empirical relation was obtained by Sho-Wei et al. [28] us¬ 
ing experimental data of longitudinal relaxation times of 
methane gas for pressures up to approximately 20 MPa. 
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FIG. 4. Longitudinal (Ti) relaxation times of bulk methane 
gas at various densities at 22 °C. The ID data were pro¬ 
cessed using the ILT algorithm [27] and each experimental 
point (blue circles) corresponds to the log mean of the entire 
ID distribution. The solid red line corresponds to the empiri¬ 
cal prediction of Ti relaxation times for spin rotation interac¬ 
tion taken from [28]. As discussed in the text, the empirical 
prediction fits the data implying that spin rotation interac¬ 
tions dominate the relaxation Hamiltonian at even such high 
densities for bulk methane. 

In Fig. 4 we compare our measured T\ relaxation times 
with the computed values using the empirical relation 
given by eq. 5 for the entire pressure range from 3.4 
MPa to 68.9 MPa. Excellent agreement (y 2 /do/ ~ 1 
[29]) was found between the empirical prediction of eq. 5 
right from lower pressures up to 68.9 Mpa (0.307 g/cm 3 ) 
demonstrating that TLtotai is dominated by the spin rota¬ 
tion interaction throughout this pressure range. It should 
be noted here that no phase transition effects were ob¬ 
served when traversing the critical pressure (P c = 4.6 
MPa) at T = 22° C as indicated by the T± and T 2 re¬ 
laxation times. Previous studies across the critical point 
(traversing both the critical pressure and temperature) 
showed changes in the relaxation times due to critical 
density fluctuations [30, 31] and distribution of correla¬ 
tion times [32], but such effects were not seen as we are 
far from the critical points in our experiments. 

B. Relaxation of methane gas in Vycor glass 

Vycor porous glass is composed of of 96% 5*02 and 4% 
Br 20 ?, and the samples used in this study had a porosity 
of 28% and an average pore size of 5.7 nm. This results 
in a high surface to volume ratio of 1.25nm -1 that leads 
to enhanced surface interaction of the confined methane 
molecules with the pore walls. The NMR relaxation of 
spin 1/2 nuclei in porous media generally occurs through 


FIG. 5. T 1 —T 2 correlation map for methane gas in Vycor glass 
at a pressure of 27.6 MPa and temperature of 22°C is shown. 
While only one peak can be resolved in the Ti dimension, the 
T 2 dimension allows for separation of different components 
in the form of three different peaks. The right most peak 
correspond to the annulus gas while the center and the left 
most correspond to the methane gas in the Vycor pore space. 
The intensity bar on the right is in [a.u]. 

two dominant pathways. The first is the slow motions at 
the surfaces resulting in enhanced relaxation caused by 
mechanisms such as Reorientation Mediated by Transla¬ 
tional Diffusion (RMTD) [33]. This effect is due to the 
geometrical characteristics of the confining system. The 
second is the dipolar relaxation due to the presence of 
paramagnetic impurities in the porous medium. As Vy¬ 
cor does not host significant quantities of paramagnetic 
impurities [34], the methane gas relaxation in this case 
is expected to be dominated by the slow motions at the 
pore walls. 

1. 2D Ti — T 2 relaxation 

The results of the T\ — T 2 correlation experiments of 
methane in Vycor at 27.6 MPa are shown in Fig. 5. The 
relaxation times were measured using an inversion recov¬ 
ery pulse sequence followed by the CPMG sequence as 
shown in Fig. 2 [35] and the results were analyzed us¬ 
ing an Inverse Laplace Transform (ILT) algorithm [22]. 
Three well distinguishable peaks were observed in the 
direct dimension (T 2 ). On the other hand, the indi¬ 
rect dimension shows a single T\ peak at a value of 
about 4.0s. This is because the T\ relaxation mecha¬ 
nism is only sensitive to motions at the Larmor frequency 
(/ = 400 MHz) and is therefore not affected by the slow 
motions of methane gas in the different environments. 

The different peaks in the X 2 dimension are a result of 
the methane molecules in different environments expe- 
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riencing different relaxation pathways. In the following 
sections we identify the origin of the different components 
of relaxation in Vycor glass and use this information to 
quantitatively assess the adsorption processes at the high 
and low pressure regimes. 

2. Spin-spin relaxation 

The pressure dependence of the T 2 relaxation times of 
methane gas in Vycor is shown in Fig. 6 . The peak 
with the longest relaxation time, labeled 7 in Fig. 6 , 
corresponds to the free methane gas in the annulus of the 
NMR tube (between the Vycor cylinder and the sample 
holder) and above the Vycor glass cylinder. This free 
methane gas relaxes mainly via spin rotation interaction 
resulting in relaxation times comparable to those of the 
bulk gas. The small difference in the relaxation times 
of these methane molecules from that of the bulk (2.4s 
instead of 2.9s at 27.6 MPa) is due to surface relaxation 
occurring on the Vycor rod surface and the NMR tube 
walls. The identity of this phase was further confirmed 
by carrying out experiments with a teflon tape covering 
the annulus space, resulting in drastic reduction of this 
long relaxing peak. 

The /? component of Fig. 6 corresponds to the methane 
gas in the pore space of Vycor. The significant reduction 
in the spin-spin relaxation time of the /3 peak with re¬ 
spect to the annulus gas (7 component) is because of the 
enhanced surface interaction of the methane molecules 
with the pore walls. The strong pressure dependent in¬ 
crease in the intensity of this peak compared to bulk 
methane gas is due to adsorption of methane molecules 
on the pore walls. Adsorbed gas has much higher density 
compared to bulk gas resulting in higher spin density in 
the pores. As the time scale of exchange between the 
free and adsorbed gas is much faster than the NMR re¬ 
laxation time scale [ 11 ], these two components cannot 
be separated and result in a single T 2 distribution for the 
free and adsorbed phases. The a component corresponds 
to the small number of methane gas molecules trapped 
in pore spaces of smaller dimensions. Vycor porous glass 
has a pore size distribution that peaks at 5.7 nrn and 
also hosts a few smaller pores due to geometrical dis¬ 
order. These molecules undergo even higher surface in¬ 
teractions resulting in relaxation times on the order of 
lms. The smaller tortuous, dead-end channels can also 
reduce molecular mobility [36]. A small background sig¬ 
nal with the same relaxation time as the a component 
was observed before gas injection which may correspond 
to trapped water molecules that could not be evacuated 
by the vacuum procedure and is shown as an inset in Fig. 
6 . 



FIG. 6. Transverse relaxation times (T 2 ) of methane gas in 
Vycor are shown for various pressures at a constant tempera¬ 
ture of 22 °C. The pressure dependence of the signal intensity 
and the T 2 values for the different components can be iden- 
tihed. The a component is zoomed in the inset to highlight 
the background signal of Vycor. 

3. 2D T 2 — I 2 relaxation 

To further understand the dynamics of the different T 2 
distributions, diffusive coupling of the methane gas be¬ 
tween the different environments was probed through the 
two dimensional T 2 — T 2 correlation measurements [37- 
39] . This experiment consists of the measurement of spin- 
spin relaxation time in the indirect dimension (“evolution 
time”) using a CPMG sequence, followed by storage of 
the magnetization in the z-axis (“storage time”), during 
which the spins relax by Tj relaxation. The latter is fol¬ 
lowed by measurement of the spin-spin relaxation time 
during acquisition using another CPMG sequence. The 
number of CPMG echoes in the indirect dimension was 
varied to acquire a 2D T 2 — T 2 plot for different storage 
times. As the spin lattice relaxation times (T) ~ 2s) at 
400 MHz are very long compared to the X 2 , the methane 
gas molecules could diffuse from one environment to an¬ 
other during the storage intervals resulting in non diag¬ 
onal peaks in the 2D T 2 — T 2 map. 

The exchange between two coupled environments may 
be quantified via the variation of the cross peak inten¬ 
sities with respect to the storage time, only in the slow 
exchange limit, i.e. the exchange is slow during the evo¬ 
lution times [37]. This condition is not satisfied in our 
experiments with significant exchange occurring between 
the coupled environments during the CPMG sequence 
(evolution time), and therefore accurate exchange times 
could not be extracted. However, the experimental data 
shown in Fig. 7 suggest that all three environments are 
coupled at the time scales of 100 ms, as cross peaks ap¬ 
pear in the T 2 — T 2 correlation plots. This indicates that 
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FIG. 7. Two dimensional map of T 2 — T 2 correlation experi¬ 
ment of methane in Vycor glass acquired with a storage time 
of 100ms. The diffusive coupling between different geomet¬ 
rical compartment results in non diagonal cross peaks. The 
first two peaks with the shorter relaxation times in the di¬ 
agonal axis correspond to the gas stored in the Vycor pore 
space while the right most peak correspond to the gas in the 
annulus space. 


the different gas species are coupled and that the a and 
/3 peaks together correspond to the gas molecules in the 
Vycor pore space and therefore should be integrated to 
obtain the total stored gas. 


C. Methane gas Hydrogen-Index in porous media 

The gas stored in Vycor glass pore space is the sum 
of the a and j3 T 2 distributions shown in Fig. 6. In Fig. 
8 the methane gas stored in the Vycor pore space as a 
function of pressure is plotted. Note that the background 
signal as shown in the inset of Fig. 6 was subtracted to 
obtain the signal from the gas stored in the Vycor pore 
space The stored gas was quantified using the calibration 
of the system to the signal from the 1 H nuclei in pure 
water. The results of the calibration of the system based 
on various volumes of pure water are shown as an inset of 
Fig. 8. It should be noted that by using a Vycor sample 
2.96 in diameter leads to a low filling factor. However, 
the results shown in Fig. 8 indicate that the low filling 
factor presents no challenge. 

The Hydrogen Index (HI) of a fluid is an important 
parameter for oil and gas exploration as it is used for 
the interpretation of various downhole logging measure¬ 
ments, including NMR, to determine the quantity of flu¬ 
ids such as hydrocarbons and water. The HI is defined as 
the density of 1 H in the sample, relative to that of water 
[8]. In the case of gas filled rock samples the hydrogen 
index of the gas phase is given by [40], 


FIG. 8. Methane gas stored in Vycor plotted as a function 
of pressure. The stored methane gas is shown to increase 
steeply at lower pressures and slowly saturates as the fluid 
density approaches that of liquid methane. The calibration 
of the signal intensity to the number of protons was carried 
out by using various volumes of H 2 O and is shown in the inset. 
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FIG. 9. Effective Hydrogen index (HI) computed based on 
the measurement of the gas stored in the Vycor pore space. 
The free methane gas hydrogen index at various pressures is 
also shown as a comparison. 


PgW'H g 

P W ■ MWg ’ 


( 6 ) 


where p w =0.11 moles- cm -3 and p g are the molar and 
mass densities of the water and gas, MW g is the molec¬ 
ular weight of methane gas and ring is the number of 
hydrogens in a single methane molecule ( njj g = 4). 

While the HI and its pressure dependence for bulk oil, 
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water and natural gas are well known [40] , the values for 
natural gas in shale nanopores are difficult to determine. 
This is because the gas in the shale pores exists as free 
and adsorbed phases whose densities also have a complex 
pressure dependence. 

For petrophysical and log interpretation applications 
an effective HI is defined which is a weighted average of 
the free and adsorbed gas hydrogen indices in the porous 
medium. The hydrogen index of bulk methane gas oc¬ 
cupying a pore volume equivalent to that of the Vycor 
sample (length of 11.50 mm, volume of 0.022 ml) is com¬ 
pared with the measured HI of methane gas in Vycor 
in Fig. 9. The effective Hydrogen Index of methane 
in the porous media is consistently higher than the 
bulk methane gas under the same conditions (pressure- 
volume-temperature) due to the increased density of the 
adsorbed gas. The measurement of the effective HI on 
real rock sample cores at formation temperatures and 
pressures using this methodology potentially enables di¬ 
rect log interpretation for obtaining fluid quantities. 

D. Monolayer and multilayer adsorption 

The gas inside the Vycor pore space exists in two 
phases, namely as adsorbed gas on the pore walls and 
as free gas in the center of the pore. As the free and 
adsorbed gas are in fast exchange in the NMR relaxation 
time scales (1ms to 50ms), the measured X 2 relaxation 
times of the gas in the pores are a weighted average of 
the two phases. The pressure dependence of the free gas 
can be approximated by that of bulk methane while that 
of the adsorbed gas is more complex. 

The quantity N excess can be defined as the excess num¬ 
ber of gas molecules in the pore space accommodated 
due to adsorption, compared to the number of bulk gas 
molecules that would have occupied the same pore vol¬ 
ume in the absence of adsorption. This quantity is given 
as [7, 41], 

-^excess = Pfree^pore? (7) 

where Nq is the total gas in the porous medium, pf ree 
is the density of the bulk methane gas and Ip 0 re is the 
volume of the pore space as specified by the manufac¬ 
turer. The total gas in the pore space corresponds to the 
sum of the a and (3 contributions of the distributions 
(Fig. 6). The second term pfreeV^re corresponds to the 
gas in the pore space assuming no adsorption. 

The excess gas in Vycor pore space may also be written 
as 

^excess = Aa ds — PfreeKds = -^ads ( 1 “ , (8) 

' Pads ' 

where V ads is the number of adsorbed gas molecules, 
Pfree and p a ds are the densities of the free and adsorbed 
gas phases and V a d s is the volume of the pore space oc¬ 
cupied by the adsorbed gas. The pressure dependence 


x 10“ 4 



FIG. 10. The quantity N excesa as a function of pressure deter¬ 
mined experimentally is fit to the modified Langmuir isotherm 
(eq. 10) at the low pressure regime. This least squares fit 
with three degrees of freedom helps determine the adsorbed 
density layer, p a d s , Langmuir pressure, Pl and the maximum 
number of adsorbed molecules, Nl- The inset figure shows 
Vex cess for the entire pressure regime, which can be described 
by multilayer adsorption models as described in the subse¬ 
quent section. 

of -/V ads helps identify the adsorption process due to its 
unique characteristics for monolayer [42] or multilayer 
adsorption [43] corresponding to the low or high pressure 
regimes. In the following sections we identify the pressure 
regimes corresponding to the Langmuir monolayer ad¬ 
sorption and BET multilayer adsorption and determine 
the density and volume occupied by the adsorbed layer, 
p ads and I4ds, as well as the number of adsorbed layers 
formed at higher pressures. 

1. Monolayer adsorption 

In the low pressure regime (0.7 Mpa< P < 39.7 MPa) 
the adsorption process is established by a monolayer 
of adsorbed molecules which can be described by the 
Lamguir isotherm given as [42], 

A^ds = Nl p + , (9) 

where P is the applied pressure, Nl the maximum num¬ 
ber of adsorbed molecules at infinite pressure and Pl the 
pressure at which half of the adsorption sites are occu¬ 
pied. 

The resulting excess number of gas molecules can then 
be described by the modified Langmuir isotherm, ob¬ 
tained by combining eq. 8 and eq. 9 as [41], 
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eq. 11 into eq. 8 we obtain the modified BET isotherm 
given as, 


N K 


= N l 


P + Pl 
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The N excess computed using eq. 7 is shown in Fig. 10 up 
to P = 39.6 MPa where the Langmuir monolayer model 
still holds. The three unknown parameters Nl, Pl and 
Pads were obtained by performing least squares fit of the 
data to eq. 10 with y 2 /do/ = 2.4 [29, 44]. The values 
obtained by the fit for the Langmuir pressure, Langmuir 
volume and the equilibrium adsorbed gas density were, 
P L = 10.80 ± 0.43 [MPa], N L = (2.16 ± 0.06) • 10- 4 [mol • 
g _1 ] and p a ds = (2.43 ±0.02) • 10 _2 [mol • ml -1 ] and agree 
well with the values reported earlier for adsorption stud¬ 
ies in activated carbon [41] . The density of the adsorbed 
layer (p a ds) was determined to be 8.5% lower than the 
density of liquid methane. 

It can be seen in Fig. 10 that the quantity N excess 
increases at low pressures to reach a maximum at ap¬ 
proximately 13.8 MPa and then decreases at higher pres¬ 
sures. This is because of the increase in the number of 
adsorbed molecules at lower pressures (P < 13.8 MPa), 
due to the presence of a large number of vacant surface 
sites. At pressures greater than 13.8 MPa the number 
of adsorbed molecules slowly plateaus out, resulting in 
a slow decrease in N excess . At higher pressures and for 
monolayer adsorption, as the density of the free gas in the 
center of the pore approaches that of the adsorbed gas, 
-^excess should asymptotically approach zero. However, 
at pressure values of P > 39.6 MPa we observe that the 
quantity N excess starts to slowly increase again as shown 
in the inset of Fig. 10. This constitutes a deviation from 
the Langmuir isotherm behavior resulting from the for¬ 
mation of multiple adsorbed layers on the pore surface. 
In the following section we analyze the entire data set up 
to 89.7 MPa using a generalized isotherm, accounting for 
the excess number of molecules in the pore space due to 
multilayer adsorption. 


2. Multilayer adsorption 

In this section we consider the Brunauer-Emmett- 
Teller (BET) model which describes multilayer adsorp¬ 
tion to describe the adsorption of methane in Vycor 
porous glass up to pressures of 89.7 MPa. The number 
of adsorbed molecules for multilayer adsorption process 
is given by [43], 


^) W+1 1 

P_\N+1 

Po> 

i-^V 

Pads J 

( 12 ) 

It should be noted that eq. 12 reduces to a form of the 
Langmuir equation in the limit of P << Pq and N = 1, 
given as, 
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N l cP 
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iV excess = N r , cJ ' ■ (1 - -BL. j, (13) 

P 0 +cP y p ads J 

which is equivalent to the Langmuir isotherm (eq. 10) 
for P L = P Q /c. 

Performing nonlinear least squares fit of eq. 12 for 
the excess gas in the pore space to the entire data set 
(shown in Fig. 11) up to pressures of 89.7 MPa we obtain, 
P 0 = 168.0T14.1 [MPa], N L = (1.4=h0.3)• 10~ 4 [mol-g” 1 ], 
c = (21.6 ± 4.6) , N = (9.7 ± 3.1). The fit had y 2 /do/ = 
0.7 [29, 44], The value for p a d s obtained for the monolayer 
Langmuir adsorption in the low pressure regime (see sec. 
IIID1) was used as a constant in eq. 13 during least 
squares minimization in order to reduce the complexity 
of the fit. 

The constant c obtained by the fit results in a (Pi — 
En) = 7.5 kJ/mol energy gap between the adsorbed lay¬ 
ers at T = 22°C, where E\ is the activation energy of 
adsorption. This value agrees with the activation energy 
values reported earlier for methane adsorption on acti¬ 
vated carbon [45]. The saturation pressure Po and con¬ 
stant c obtained by fitting the experimental data help 
determine the Langmuir pressure (Pl = Po/c). The fit 
value for the Langmuir pressure is in good agreement 
with that obtained by the monolayer Langmuir model 
for the low pressure regime (sec. Ill D 1-Monolayer ad¬ 
sorption). The number of layers of methane with average 
molecule diameter of 3.9 A that can be formed into the 
Vycor pore with an average diameter of 5.7 nm is approx¬ 
imately 7.4. Note that the fit overestimates the number 
of adsorbed layers (N = (9.7 ± 3.1)) formed in the Vy¬ 
cor pore but the obtained theoretical value is within the 
experimental uncertainty. 


Ar AT cx (1 — (AT + l)x N + Nx n+1 ^ 
ads = Nl (T^±) V 1 + (c - 1 )® - cx N +' )' 


( 11 ) 


3. Separation of free and adsorbed gas 


In the above equation, x = P/Pq, where Po is the 
saturation pressure, N is the number of adsorbed lay¬ 
ers, Nl is the number of molecules in a monolayer and 
c oc e ( E i~ E N)/RT j g a un ttless constant related to the heat 
energies of the first and N layers (E\,En). Substituting 


To separate the gas in the Vycor nanopores into free 
and adsorbed fractions for the entire pressure range from 
0.7 MPa-89.7 MPa we make use of the density of the ad¬ 
sorbed layer obtained from the Langmuir isotherm fit. At 
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FIG. 11. The quantity lV ex cess is fit for the entire pressure 
range using the modified multilayer BET adsorption model 
(eq. 12). This least squares fit with four degrees of free¬ 
dom enables the determination of the adsorbed layers N, the 
constant c, the saturation pressure Pq and the number of ad¬ 
sorbed molecules in a single layer Nl . The value for p a d a 
obtained from the monolayer Langmuir adsorption fit in the 
low pressure regime (fig 10) was used as an input in the fit. 

low pressures (P < 39.6 MPa) the N eKcess can be alter¬ 
natively written in terms of the density of the adsorbed 
layer as, 

^excess |p=39.6MPa = Pads Fids Pbulk |p=39.6MPa Fads- 

(14) 

From the above equation the volume of a monolayer 
of adsorbed gas (V^ds) can be determined using the p a d s 
obtained from the Langmuir fit and the known values 
°f Pbuik- Using the adsorbed gas volume together with 
the pore volume of Vycor, the gas stored in the Vycor 
nanopores can be separated into free and adsorbed gas 
on the pore walls, for pressures up to 39.6 MPa. For pres¬ 
sures higher than 39.6 MPa, V a< i s is pressure dependent as 
multiple layers start to form and is therefore determined 
as, 

at ( p \ |P=89.7MPa 

( TD\ |P=89.7MPa_ excess V/ ) lp=39.6MPa ( 
^adslp; IP=39.6MPa * l lc V 

Pads Pbulk 

The results of the free, adsorbed and total methane 
gas in Vycor porous glass are shown in Fig 12. Up to 
a pressure of 39.6 MPa the “free gas” molecules in the 
pore interiors (red circles in Fig 12) show a bulk-like pres¬ 
sure dependence while the adsorbed gas fraction increases 
faster at lower pressures but slowly saturates at higher 
pressures. This indicates that at low pressures (P < 12 
MPa) more gas is stored in the pore space as adsorbed 
gas. As discussed in Sec. Ill D 2, multiple adsorbed layers 
start to form for pressures P > 39.6 MPa and therefore 



FIG. 12. Separation of the total gas into free and adsorbed. 
The blue squares correspond to 1 H inside Vycor. Based on 
the I4da calculated from eq. 14, the total number of adsorbed 
1 H was calculated (black diamonds). The red circles corre¬ 
spond to the total number of free 1 H in the pores. The bulk 
methane gas in Vycor (assuming no adsorption) is also shown 
in magenta triangles. 

a rapid drop is observed in the “free gas” fraction. On 
the other hand the adsorbed gas in the pore space shows 
a rapid increase up to P = 89.7 MPa, where it is approx¬ 
imately 91% of the total methane gas stored in the pore 
space. The amount of bulk methane gas that could be 
accommodated in the pore volume, in the absence of ad¬ 
sorption is also shown in magenta. This way the total gas 
stored in the Vycor nanopores can be separated into free 
and adsorbed fractions in the entire pressure range from 
0.7 Mpa to 89.7 MPa using high field NMR relaxometry. 

IV. CONCLUSIONS 

In the present work we demonstrate the application of 
high field NMR for characterization of the storage prop¬ 
erties of methane gas in nanometer sized pore spaces. 
We carry out laboratory experiments at a wide pressure 
range (0.7 MPa to 89.7 MPa), that covers most condi¬ 
tions encountered in naturally existing shale gas and tight 
gas reservoirs. We demonstrate that the relaxation of the 
bulk methane gas, which forms the dominant fraction of 
natural gas, is dominated by the spin rotation Hamilto¬ 
nian up to the highest studied pressures of 89.7 MPa. 
The sensitivity of the X 2 relaxation distributions to the 
slow motions enable the separation of the gas molecules 
in the pore space and thereby aid gas storage and HI 
measurements as a function of pressure. 

The quantity of methane gas in Vycor pore space was 
used to determine the excess gas due to adsorption on 
the pore walls. Two pressure regimes that correspond 
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to monolayer and multilayer adsorption were identified. 
In the low pressure regime from 0.7 < P < 39.6 MPa 
the excess gas was fit to the Langmuir isotherm and the 
Langmuir pressure, volume and equilibrium density of 
the adsorbed layer were determined. The adsorbed gas 
density was discovered to be 8.5% lower than that of liq¬ 
uid methane. In the high pressure regime from 39.6 Mpa 
to 89.7 MPa the pressure dependence of the excess gas 
shows evidence of multilayer adsorption and was mod¬ 
eled using BET theory to determine the number of ad¬ 
sorbed layers formed. Using the adsorbed gas density, 
the gas in the Vycor pore space was separated into free 
and adsorbed fractions for the entire pressure range. The 
adsorbed gas fraction was shown to increase sharply for 
pressures above 39.6 MPa due the formation of multiple 
adsorbed layers, while the quantity of free gas molecules 
in the pore space decreased correspondingly. 

Vycor porous glass with 5.7 nm pore size used in 
these studies can be considered as a model system 
to understand storage properties of natural gas in 
shale gas and tight gas rock formations. While the 
pore size of the Vycor porous glass is comparable to 
the small pores in the shale rock samples, important 
differences include the low porosities, permeability and 


also the presence of paramagnetic impurities in the 
naturally occurring shale rocks. The presence of internal 
gradients due to paramagnetic impurities in the shale 
nanopores are averaged out by the fast diffusing gas 
molecules, but can affect the annulus gas outside the 
rock resulting in their overlap in the relaxation distri¬ 
butions. The low porosity and permeability of shale 
rocks also make the NMR experiments time consuming 
in order to achieve acceptable signal to noise ratio and 
complete gas saturation. In spite of these challenges, 
the sensitivity and precision of the high field NMR 
relaxometry experiments make it a viable methodology 
for studying the gas storage and transport character¬ 
istics of natural gas in such micro and mesoporous media. 
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